Abstract 



Self-adjoint Schrodinger and Dirac operators with Aharonov-Bohm 

and magnetic-solenoid fields 

o 

O ! D.M. Gitman* A. Smirnov, I.V. Tyutin] and B.L. Voronov* 

> 
O 

We study all the s.a. Schrodinger and Dirac operators (Hamiltonians) both with 
pure AB field and with magnetic-solenoid field. Then, we perform a complete spectral 
analysis for these operators, which includes finding spectra and spectral decomposi- 
tions, or inversion formulas. In constructing the Hamiltonians and performing their 
spectral analysis, we respectively follow the von Neumann theory of s.a. extensions 
^ | of symmetric differential operators and the Krein method of guiding functionals. The 

examples of similar consideration are given by us in arXiv:0903.5277, where a non- 
relativistic particle in the Calogero potential field is considered and in Theor. Math. 
Phys. 150(1) (2007) 34, where a Dirac particle in the Coulomb field of arbitrary charge 
is considered. However, due to peculiarities of the three-dimensional problems under 
consideration, we elaborated a generalization of the approach used in the study of the 



■ Dirac particle. 

7-H ■ 

o\ '. 1 Introduction 

O 

Aharonov-Bohm (AB) effect [I] plays an important role in quantum theory refining the status 
^ ■ of electromagnetic potentials in this theory. First this effect was discussed in relation to a 

study of interaction between a non-relativistic charged particle and an infinitely long and 
infinitesimally thin magnetic solenoid field (further AB field) which yields a magnetic flux 
<3> (a similar effect was discussed earlier by Ehrenberg and Siday [2]). It was discovered that 
particle wave functions vanish at the solenoid line. In spite of the fact that the magnetic 
field vanishes out of the solenoid, the phase shift in the wave functions is proportional to 
the corresponding magnetic flux [3j. A non-trivial particle scattering by the solenoid is 
interpreted as a possibility for quantum particles to "feel" potentials of the corresponding 
electromagnetic field. Indeed, potentials of AB field do not vanish out of the solenoid. For 
the first time, a construction of self-adjoint (s.a. in what follows) Schrodinger operators with 
the AB field was given in [I]. First, the need for s.a. extensions of the Dirac Hamiltonian 
with the AB field in 2 + 1 dimensions was recognized in El [7j. S.a. extensions of the Dirac 
Hamiltonian with the AB field in 3 + 1 dimensions were found in [8] , see also [9j \TU\ . The 
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physically motivated boundary conditions for the particle scattering by the AB field and 
a Coulomb center were studied in [TT1 [T2] . A splitting of Landau levels in a superposition 
of parallel uniform magnetic field and AB field (further magnetic- solenoid field) gives an 
example of AB effect for bound states. First, exact solutions of Schrodinger equation with 
the magnetic-solenoid field (non-relativistic case) were studied in [13]. Exact solutions of 
the relativistic wave equations (Klein-Gorgon and Dirac) with the magnetic-solenoid field 
were obtained in [TH [151 05] and used then to study AB effect in cyclotron and synchrotron 
radiations, see [TBI fT6l fT7] . Later on the problem of self-adjointness of the Dirac Hamiltonian 
with magnetic-solenoid field was studied in [TH [19] . 

In the present work, we construct systematically all the s.a. Schrodinger and Dirac 
operators both with pure AB field and with magnetic-solenoid field. Then, we perform a 
complete spectral analysis for these Hamiltonians, which includes finding spectra and spectral 
decompositions, or inversion formulas. In constructing the Hamiltonians and performing 
their spectral analysis, we respectively follow the theory of s.a. extensions of symmetric 
differential operators [201 |2TJ, [25] an d the Krein method of guiding functionals [201 |2T]. The 
examples of similar consideration are given in [22] where a nonrelativistic particle in the 
Calogero potential field is considered and in [23] where a Dirac particle in the Coulomb field 
of arbitrary charge is considered. However, due to peculiarities of the three-dimensional 
problem under consideration, we use a necessary generalization of the approach [23]. 

We recall that the AB field of infinitely thin solenoid (with the constant flux <3>) along 
the axis z = x 3 can be described by electromagnetic potentials A AB , p, = 0, 1, 2, 3, 

^ab = (0, A AB ) , A AB = (A k AB , k = 1, 2, 3) , A 3 AB = 0, 
! $siny 2 $cosy? 

Aab ~ Aab ~ ~2^T' 

where p, if are cylindrical coordinates, x 1 = p cos if, x 2 = psimp. The magnetic field of AB 
solenoid has the form Bab = (0, 0, -Bab)- It is easy to see that outside the z axis the magnetic 
field Bab = rotAAB is equal to zero. Nevertheless, for any surface S with a boundary L being 
any contour (even an infinitely small one) around the z axis, the circulation of the vector 
potential along L does not vanish and reads § L Ap&dl If one interprets this circulation 
as the flux of the magnetic field Bab through the surface E, 

B AB da = f A AB d\ = $, 

then we obtain an expression for the magnetic field, 

Bab = ^(x^ix 2 ) , 

where the term infinitely thin solenoid comes from. 
In the cylindrical coordinates, we have 

^AB = -0P _1 Sin (f , ^AB = COS if , (f) = $/$ , 

where $o is a fundamental unit of magnetic flux, 



$ = 2-Kch/e = 4, 135 x 10 7 Gauss ■ cm 2 
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(we recall that e > is the absolute value of the electron charge). 

The magnetic-solenoid field is a superposition of a constant uniform magnetic field of 
strength B directed along the axis z and the AB field with the flux $ in the same direction. 
The magnetic-solenoid field is given by electromagnetic potentials by potentials = (0, A), 
A = (A k , k — 1,2, 3) of the form 

A 1 = A\ B -^A 2 = A\ B + ^, A* = U. (1) 

The potentials ([I]) define the magnetic field B of the form 

B = (0,0, B + B AB ). 

In the cylindrical coordinates, the potentials of the magnetic-solenoid field have the form 

—A 1 = -0p _1 sin ip , —A 2 = 0p _1 cos <p , A 3 = 0, 
cn cn 

= 0+^, 7 = ^>0, e B = sign£. (2) 
For further consideration, it is convenient to introduce the following representation: 

4> = £b (0o + A*) j 0o = Kb0] e Z, /i = e B 4> - O , < /x < 1. (3) 

The quantity /i is called the mantissa of the magnetic flux and, in fact, determines all the 
physical effects in the AB field, see e.g. 



2 S.a. Schrodinger Hamiltonians 

In this section, we consider two-dimensional and three-dimensional nonrelativistic motions 
of a particle of mass m e and charge q = e q e, e q = signg = ±1 ( positron or electron) in 
the magnetic-solenoid field. The canonical formulation of the problem is the following. The 
starting point is the "formal Schrodinger Hamiltonian H" with the magnetic-solenoid field 
that is respectively a two- or three-dimensional s.a. differential operation well-known from 
physical textbooks. In three dimensions, it is given by 

6 = 2^T ( p ~ l A ) 2 ' p = ~* W ' V = {dx ' 9yj dz) ' (4) 

It is convenient to represent if as a a sum of two terms, H 1 - and if", 

H = H 1 - + if 11 , 

where the two-dimensional s.a. differential operation if- 1 , the "formal two-dimensional 
Schrodinger Hamiltonian" with the magnetic-solenoid field, 



H 1 = M" 1 ^, 7f ± = ( -iV ± - 4 A 

V cn 




M = 2m e h~\ V x = (d x , d y ) , A ± = (A\ A 2 ) , (5) 
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A 1 and A 2 are given by (121) , corresponds to a two-dimensional motion in the xy plane 
perpendicular to the z axis, while the one-dimensional differential operation w\ 

~2 

#n = # = ^*_ p z = -ihd z , 

corresponds to a one-dimensional free motion along the z-axis. 

The problem to be solved is to construct s.a. nonrelativistic two- and three-dimensional 
Hamiltonians H and H associated with the respective s.a. differential operations H 1 - and 
H and to perform a complete spectral analysis for these operators. 

We begin with the two-dimensional problem. We successively consider the case of pure 
AB field, with B = 0, and then the case of the magnetic-solenoid field, with B ^ In the 
subsequent subsection, we generalize obtained results to three-dimensions. 



2.1 Two-dimensional case 
2.1.1 Reduction to radial problem 

In the case of two dimensions, the space of particle quantum states is the Hilbert space 
ft = L 2 (R 2 ) of square-integrable functions 4>{p), p = (x, y), with the scalar product 



^2) = / iPi(p)*p2(p) d P, dp = dxdy = pdpdip. 

A quantum Hamiltonian H 1 - should be defined operator in this Hilbert space. It 

is more convenient to deal with a s.a. operator Ti 1 - = M H 1 - associated with s.a. differential 
operation H 1 - = M H 1 - defined in ([5j). 

The construction is essentially based on the requirement of rotation symmetry which 
certainly holds in a classical description of the system. This requirement is formulated as 
the requirement of the invariance of a s.a. Hamiltonian under rotations around the solenoid 
line, the z axis. As in classical mechanics, the rotation symmetry allows separating the polar 
coordinates p and (p and reducing the two-dimensional problem to a one-dimensional radial 
problem. 

The group of rotations 50(2) in IR 2 naturally acts in the Hilbert space Sj by unitary 
operators: if S G SO (2), then the corresponding operator Ug is defined by the relation 

(U s iP)(p) = HS- 1 p)^ef). 

The Hilbert space is a direct orthogonal sum of subspaces S) m , that are the eigenspaces 
of the representation Us, 

where 9 is the rotation angle corresponding to S. 

It is convenient to change indexing, m — > I, S) m — >• Sji, as follows 

m = e (0 O — I) , I = 4> — em, 



where 



e = e q e B , (j) = [e B (j>}- 



(6) 



5 



We define a rotationally-invariant initial symmetric operator 7i ± associated with the 
differential operation Ti, 1 - in the Hilbert space = L 2 (R 2 ) as follows: 



H- 



D n x = {ip{p) : ip e ^(M 2 \{0})> 



where £>(IR 2 \{0}) is the space of smooth and compactly supported functions vanishing in a 
neighborhood of the point p =0. The domain D n ± is dense in Sj and the symmetricity of 
H 1 - is obvious. 

In the polar coordinates p and ip, Ti, 1 - becomes 



H 1 = -d 2 - p-% + p- 2 {id v + e,0) 2 , (7) 



where <f> is given by (J2J). 
For every Z, the relation 



(Sif)(fiM = ^^*°- lyp fi{p) (8) 



determines a unitary operator S 1 ; : L 2 (R + ) — ► where L 2 (R + ) is the Hilbert space of 
square-integrable functions on the semi-axis R + with scalar product 



(/,<?)= / f{p)g{p)d P . 

For every Z, we define the linear operator V\ from ij to L 2 (R + ) by setting 

(Vrt)(p) = r^tfe-****^"**. (9) 



'2ti 



o 



If ip = J^iez^i e ^> then we have ^ = SjVj^ for all I. In other words, VJ = SfPi, where Pi 
is the orthogonal projector onto the subspace S)i- However, we prefer to work with V\ rather 
than P\ because the latter cannot be reasonably defined in the three-dimensional case, where 
the Hilbert state space should be decomposed into a direct integral instead of a direct sum 
(see below). Clearly, V t if) E V(R+) for any ip e £>(R 2 \ {0}), and it follows from © and 
that 

V{H^if) = h(l)Vfll>, ip E V(R 2 \ {0}), (10) 
where the symmetric operators h(l) in L 2 (R + ) defined on D^m = T>(R + ), where it acts as 



h{l) = -d 2 p + p~ 2 (/ + /i + 7p 2 /2) 2 -l/4 



'11} 



with p = €b4> — 0o defined in fl3]). 

In view of ffTOl . for any if> G X>(R 2 \ {0}), the .^-component (H ± if))i °f 7~L~ L 4' can De 
written as 

(H^ = SiViH^ip = Sih^S^SiViif) = Sih^S^. (12) 
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Suppose we have a (not necessarily closed) operator / in Sji for each /. We define the 
operator 

in fj by setting 

H = ^2 Mi, V> = J^i- 

lei lez 

The domain Df of / consists of all ip = J^iez^i e ^ such that ip t <E Df t for all / and the 
series J^ ZgZ fiipi converges in f). The operator / is closed (self-adjoint) if and only if all fi 
are closed (resp., self-adjoint). For every I, we have Df t — Df fl fii. 

We say that a closed operator / in fj is rotationally invariant if it can be represented in 
form ffl3l) for some family of operators fi in Sji. 

By ( TT2|) . the direct sum of the operators Sih^Sf 1 is an extension of 

n 1 - c J^Sj^ost 1 . ( 14 ) 

Let ^ e (0 De s - a - extensions of the symmetric operators h(l). Then the operators 

H±(l) = SMQSr 1 (15) 
are s.a. extensions of Sih(l)Sf , and it follows from ffl4|) that the orthogonal direct sum 

is a rotationally invariant s.a. extension of the initial operator Ti, 1 . 

Conversely, let be a rotationally invariant s.a. extension of TC . Then it has the 
form (fl6|) . where 7~C^(l) are s.a. operators in Sji. Set /i e (Z) = S^H^r^Si. For all /, h e (l) are 
s.a. operators in L 2 (R+). If f e V{+), then ft/ E V(R 2 \ {0}) nJ5i and flU} and (USD imply 
that 

S$(0/ = SMQS^SJ = H^Srf = H^SJ = H^IW = SMl)f. 

Hence, h(l)f = h t (l)f, i.e., h t (l) is a s.a. extension of h(l). We thus conclude that Ti^ can 
be represented in form (fT6l) . where "H^il) are given by ffl5|) and /i c (Z) are a s.a. extensions 
of h(l). 

The problem of constructing a rotationally invariant s.a. Hamiltonian TCj- is thus reduced 
to constructing s.a. radial Hamiltonians h e (l). 

We first consider the case of pure AB field where B = 0. 

2.1.2 S.a. Hamiltonians with AB field 

In this case, we have 7 = 0, and s.a. radial differential operation h(l) (1111) becomes 

h(l) = -d 2 p + ap~ 2 , a = xf-l/4, x x = \l + I 6 Z. (17) 

It is easy to see that this differential operation and the corresponding initial symmetric 
operator h (I) are actually identical to the respective operation and operator encountered in 
studying the Calogero problem, see [22]. We can therefore directly carry over the previously 
obtained results to s.a. extensions of h (I). 
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First region: a > 3/4, In this region, we have (Z + p) 2 > 1, which is equivalent to 

Z > 1 — p or Z < — 1 — /x. 

Because / G Z and < p < 1, we has to distinguish the cases of /x = and /x > 0: 

^ = : I < -1 or I > 1, i.e., Z + 0, 
/x > : I < -2 or I > 1, i.e., Z ^ 0, -1. 

For such Z, the initial symmetric operator h (Z) has zero deficiency indices, is essentially 
s.a., and its unique s.a. extension is Zx e (Z) = h(i) (Z) = h + (Z) with the domain 

Dfan = {V>* : ip[ are absolutely continuous (a.c.) in R + , ^ Zx (Z) ^ G L 2 (M + )}. 

The spectrum of Ztm (Z) is simple and continuous and coincides with the positive semiaxis, 
spec/i(i) (Z) = M + . 

The generalized eigenfunctions Us, 

U £ (p) = (p/2) 1/2 J >tl (VSp), h {1) {l)U £ = £U £ , SeR + , 
of Zx(i) (Z) form a complete orthonormalized set in each Hilbert space Hi. 

Second region: —1/4 < a < 3/4 In this region, we have < (Z + p) 2 < 1, which is 
equivalent to 

— xx<Z<l — /xor — 1 — /x < Z < —/x . (18) 

If /x = 0, inequalities (ITS]) have no solutions for I G Z. If /x > 0, these inequalities have 
two solutions Z = Z a , where, for brevity, we introduce the notation 

l a = a, a = 0, —1. 

So, in the second region, we remain with the case of /x > 0.. 

For each I = l a (a = 0,-1) there exists a one-parameter [/ (l)-family of s.a. Hamil- 
tonians h t (l a ) = h\ a (l a ) parametrized by the real parameter A a G S (— ir/2, vr/2), where 
S (— 7r/2, 7r/2) denotes an interval with identified ends (a circle). These Hamiltonians are 
specified by the asymptotic s.a. boundary conditions at the origin, 

^ Aa (p) = c [(M 1/2+ " cosA a + (koP) 1/2 ~^ sin Aj + 0(p 3 / 2 ), p - 0, (19) 
Dh M = i> G Zfyy (K+) , V satisfies (20) 

where x a = x; a = |/x + a|, < x a < 1, and c is an arbitrary constant, whereas ko is a 
constant of dimension of inverse length. 

For A a G" (— 7r/2,0), the spectrum of each of h\ a {l a ) is simple and continuous and 
specZz Aa (l a ) = E + . 
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The generalized eigenfunctions Us, 



Ue(p) 



P 

2Q a 



Q a = l + 2A a {£/A)"" cos(7rx a ) + ( A a ) (£/4) 2 ^ > 0, 

A a = r(l - x a )r -1 (l + x a ) tanA a 
hAQU £ = £U £ , £gM + , 



(21) 



of the Hamiltonian h\ a (l a ) form a complete orthonormalized set in each Hilbert space L\ a . 

For A a G (—7r/2,0), the spectrum of each of h\ a (l a ) is simple, but in addition to the 
continuous part of the spectrum, there exists one negative level = — 4/cQ|A a |~^ a , such 

that spec/z Aa (l a ) = K+ U {e^}. 

The generalized eigenfunctions Us of the continuous spectrum, £ > 0, are given by the 
same (I2TI) . while the eigenfunction U^~' corresponding to the discrete level is 



\ 



2p 



4~ ] 



sm(7rx a ] 



7TX„ 



they together form a complete orthonormalized set in each Hilbert space Ci a . 

Third region: a = — 1/4 In this region, we have I + fi — 0. If p, = 0, this equation has 
a unique solution I = l = 0, while if p > 0, there are no solutions, and we remain with the 
only case of p, = 0. 

For I = Iq, there exists a one-parameter £/ (l)-family of s.a. Hamiltonians h e (Iq) = 
h\ {Iq). parametrized by the real parameter A G § (— 7r/2, 7r/2)) These Hamiltonians are spec- 
ified by the asymptotic s.a. boundary conditions at the origin 

^ A (p) = c [p 1/2 In (kqp) cos A + p 1/2 sin A] + 0(p 3/2 hip), p -> 0, (22) 

( the constants c and ko are of the same meaning as in ([19]) ). and their domains are given by 

D hx{ i o) = : V e ^ (lo) (»+) , V> satisfies (J22])}, (23) 



(R+) = : V*, V>* are a.c. in E + , V*, & (fo) V, e £ 2 (» + )}. 



The spectrum of /ia (lo) is simple. For |A| = ir/2, the spectrum is continuous and non- 
negative, spech± n /2 (lo) = M+. For |A| < ir/2, in addition to the continuous part of the 
spectrum, £ > 0, there exists one negative level £[ ^ = — 4kq exp [2 (tan A — C)], where C is 
the Euler constant, such that 



spechx (lo) = \ £ 



•(-) 



U 



The generalized and normalized eigenfunctions Us of the continuous spectrum are 

AJ (ye P ) + ^n (yep 



Ue(p) 




p 

2 ( A 2 + 7r 2 /4 



A = tan A - C - In (V£/2k ) , h x (Z ) = , £ G R+, |A| < tt/2, 
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while the normalized eigenfunction ) corresponding to the discrete level is 

U^-\p) = ^2p\£^\K (y^Flp) , h (Z ) t/ ( - } = ^[/^ , |A| < tt/2, 
they together form a complete orthonormalized set in the Hilbert space £/ . 

Complete spectrum and inversion formulas In the previous subsubsecs., we con- 
structed all s.a. radial Hamiltonians h e (I) associated with the s.a. differential operation 
h (I) as s.a. extensions of the symmetric operator h (I) for any Z G Z and for any any <p 
and \i. We assemble our previous results into two groups. 
For \i = 0, we have 

h t (l) = h (1) (l), Z^Zo, D fc(1)(0 = (R+) , 

k{lo)=h x {lo) A G § (-tt/2, tt/2) , (24) 

the domain D hx ^ ) is given by (123]) ; 
For /x > 0, we have 

Zi c (Z) = (Z) , Z ^ Z a = a = 0, -1, D h(1)(J) = Dl (l) (R + ) , 

(l a ) = h a (L) , A Q G § (-tt/2, tt/2) , (25) 

the domain D hx is given by (|20|) . 

As a final result, we find a family of all s.a. rotationally- invariant two-dimensional nonrel- 
ativistic Hamiltonians = M _1 7i^ associated with the s.a. differential operation H 1 - (jHJ) 
with B = 0. Each set of possible s.a. radial Hamiltonians h t (I) generates a s.a. rotationally- 
invariant Hamiltonian in accordance with the relations f|T5|) and ffT6l) . 

When presenting the spectrum and inversion formulas for H^~, we also consider the case 
of n = and the case of > separately . We let E denote the spectrum points of 
and let ^ e denote the corresponding (generalized) eigenf unctions. The spectrum points of 
the operators h e (l) and are evidently related by £ = ME. Therefore, when writing 
formulas for eigenf unctions ^>e of the operator in terms of eigenf unctions Ug of the 
operators h c (l), we has to introduce the factor \j ^/2~np e ieq ^°~ l ' v in accordance with eq.(JH]) 
with e = e q (because e# = 1), to make the substitutions 8 = ME and £^~' = ME^J, 

^ = ME X ^ for the respective points of the continuous spectrum and discrete spectrum, 
and, in addition, to multiply eigenfunctions of the continuous spectrum of the operators h t (I) 
by the factor \[M because of the change of the spectral measure dS to the corresponding 
spectral measur^ dE. 

For /i = 0, there is a family of s.a. two-dimensional nonrelativistic Hamiltonians 
parametrized by the real parameter A G § (—tt/2, tt/2) , Hj- = 

^ = YT^(l)@Hi(k), 

H 1 - (l) = M- 1 S l h {1) (l)Sf\l^l , 
Hi(l ) = M- 1 S lo h x (l )S h l . 

^rom the physical standpoint, the latter is related to the change of the "normalization of the eigenfunc- 
tions of the continuous spectrum to S function" from 8(£ — £') to S(E — E'). 
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The spectrum of Hj~ is given by 



specify 



E ( ~ ] = -4M- 1 ^ eX p [2(tanA - C)] , |A| < tt/2 
0, A = ±tt/2 



u: 



The complete set of orthonormalized eigenfunctions of consists of the generalized 
eigenfunctions ^i >E (p) of the continuous spectrum, E > 0, 

%, E (p) = (M/47r) 1/2 e^°-0^ (VMEpj , I ^ l , 



M 



^ 471 [\ + 7T 2 /4^ 

A = tan A - C - In (^/ME/2K^j , 



AJ r 



+ -n (Vme p ) 



and the eigenfunction ^(p) corresponding to the discrete level e[ ^ in the case of |A| < tt/2, 



(-) 



such that 

= E%,e(p), H^lM = E*l E (p), E>0, 
The corresponding inversion formulas are 



*(p)= £ / * l (E)* ltE (p)dE+ $ lo (E)*l E (p)dE + <S> lo ^ o (p), 



HE) = J dp%, E (p)*(p), $ i0 (S) = | dp^ oE (p)*(p), * Io = J dp**(p)*(p), 
dp|^(p)| 2 = E / |$,(^)| 2 ^+l^o| 2 > V^GL 2 



the terms including and ^(p) are absent in the case of |A| = tt/2. 

For fi > 0, there is a family of s.a. two-dimensional nonrelativistic Hamiltonians Hj 

r-L 



parametrized by two real parameters A a G 8 (—tt/2, tt/2)) , if^ = Hh a = 0, —1, 



lGZ,l^la a 

H L (l) = M- 1 S l h {1) (l)Sf 1 ,l^l a , 
Hi(la) = M- 1 S L h Xa (l a )Si a \ 
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The spectrum of Hf\, is 
specH(- Xa} = 



Eg = -AM-^Xat*" 1 , \ a e (-vr/2,0) 
0, A a ^(-7r/2,0) 



U 



where x a = \p + a\ , A a = T(l - x a )T x (l + x a ) tanA a . 

The complete set of orthonormalized eigenfunctions of H^~ x } consists of the generalized 

eigenfunctions ^/^(p), I ^ l a , and ^° E (p) of the continuous spectrum, E > 0, 



*, )B (p) = (M/4tt) 1/2 e ^<HV [4ME^ ,x l = \l + p\,l^l a , 

J« a (yfMEp) + X a (yME/2K y >Ca J_„ a (y/MEp 



Je q ((f> -l a )<p 



4vrQ t 



Q a = 1 + 2\ a {ME /4)** cos(7rx a ) + A a (ME/4) 



and the eigenfunctions ^ a (p) corresponding to the discrete levels ■* in the case of A a G 
(-tt/2,0) 



such that 





M 2 




sin(7rx a ) 




7T 2 




H x % E (p) = 




I? la, 




= e[~ 


"Wo 

l b 


[P), b = 



ME 



(") 



The corresponding inversion formulas are 

/oo r /*oo 

^(E)^ l>E (p)dE + J2 / ^ la (E)^t E (p)dE + ^X:(p) 

$,(£7) = / dp%M^f{ P ), Z + l a , 

$ la (E) = J d P vfc E (pm P ), $ la = J d P *f;(p)*(p), 

//*oo 
dpi^(p)i 2 = ^ / i$ z (^)i 2 ^+Xli^i 2 ' WgL 



- 2 /Tn>2\ 



the terms with $/ a and ^"(p) are absent in the case of A a ^ (—tt/2,0). 
We now consider the case of the magnetic-solenoid field where B ^ 0. 

2.1.3 S.a. Hamiltonians with magnetic-solenoid field 

In this case, the radial differential operation h (I) is given by (iTTj) with 7 = ^IM ^ 0, or 



h(l) = -3f; + g 1 p-* + g 2 p* + £V>, 



9l = ni - 1/4, x, = \l + p\ , g 2 = 774, = l(l + p) 
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Up to the constant term 8\ , this s.a. differential operation is identical to the one- 
dimensional Schrodinger operation —d 2 x + g\X~ 2 + g 2 x 2 studied by us recently (to be pub- 
lished). We therefore directly carry over the obtained results to s.a. extensions of h(l). 
We note that as in the case of pure AB field, a division to different regions of giis actually 
determined by the same term g\p~ 2 = ( x 2 — 1/4)/ p~ 2 singular at the origin and independent 
of the value of B. 



First region: g\ > 3/4 In this region, we have (I + /i) 2 > 1 and as before, we distinguish 
the cases of p = and /i > : 

fM — : I < — 1 or I > 1, i.e., I ^ l , 
/i > : I < -2 or I > 1, i.e., / ^ l a . 

For such /, the initial symmetric operator h (/) has zero deficiency indices, is essentially 
s.a., and its unique s.a. extension is h c (I) = (/) = h + (/) with the domain D* h ^ . 

The spectrum of h(i) {I) is simple and discrete, 

specfyi) (I) = {£ l>m = 7 (1 + \l + fi\ + (I + /I) + 2m) , m e Z + } . (26) 
Eigenf unctions of the Hamiltonian (I) are 

U&ip) = Ql,m (7/2) 1/4+ " /2 p 1/2+ "e^ 2 / 4 $(-m, 1 + x i;7P 2 /2), 
^tT(1 + x ; + m)^ 1/2 



Q 



Lm 



m\T 2 (l + hi) 

ki)Q)u£i = £i, m u$i> ™ eZ +> ( 27 ) 

they form a complete orthonormalized set in each Hilbert space 

Second region :— l/4<a<3/4 In this region, we have < (I + p) 2 < l,or equivalently 
([TBI . We know that if fi = 0, these inequalities have no solutions for I e Z, while if \i > 
there are the two solutions I = l a = a, a = 0, —1. Therefore, we again remain with the case 
of yU > 0. 

For each I = l a , there exists a one-parameter U (l)-family of s.a. Hamiltonians h e (l a ) = 
h\ a (l a ) parametrized by the real parameter A a G § (— tt/2, 7r/2) These Hamiltonians are 
specified by the asymptotic s.a. boundary conditions at the origin 



^a/t72/°) sin + (VtT^p) cosA c 



+ 0(p 3 / 2 ), (28) 



where x a = + a|, < x a < 1, and c is an arbitrary constant! 2 ], and their domains are given 
by 

D hXa (i a ) = G D^ (Jb) (R + ) , V> satisfies ((28])}, (29) 
D^ (Jb) (R+) = : V*, are a.c. in R+, V*, (k) ^* £ L 2 (R + )}. 



2 In comparison with (119p . we fix the dimensional parameter ko by fco = v/7/2. 
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The spectrum of h\ a (l a ) is simple and discrete and is bounded from below, 

Spech Xa (L) = |^a,m = Ta,m + £^\ TU G Z+ j , 

where r ajm are solutions of the equation uj\ a (r a ^ m ) = 0, 

u\ a (W) =u+(W) sinA a + u;_(W) cosA a , 

u±(W) = r(l ± x tt )/r(l/2 ± x a /2 - W/2 7 ). (30) 
The eigenfunctions U\ a m of the Hamiltonian /iA a (/ a ) are 

^fffJ/ 3 ) = Q<hm l U +{P'i T a,m) Sm A a + M_(p; T a , m ) COS A a ] , 

1/2 



Q a , m = * a Y7 , , u Xa (W) =u + (W)cos\ a -uj4W) S mX a , 



u±(p; W) = ( 7 /2) 1/4± ^ /2 p 1 / 2± -e-^ 2 / 4 $(l/2 ± x a /2 - W/2-y, 1 ± x a ; 7 p 2 /2), 

K (Q U£] m = Sa,mU® m , (31) 

they form a complete orthonormalized set in each Hilbert space L\ a . 

Explicit expressions for the spectrum and eigenfunctions can be easily obtained in the 
cases of A a = ±7r/2 and A a = 0. In the case of A a = ±7r/2, they are given by the respective 
formulas (|2"oT) and (12TI) with the substitutions I — > l a and — > x a . In the case of A a = 0, 
these formulas are modified by the additional substitution x a — > — x a . 

Third region: a — — 1/4 In this region, we have / + p = 0. As we know, we remain with 
the only case of p = and with I — l — 0. 

For I = Zo, there exists a one-parameter U (l)-family of s.a. Hamiltonians h e (1$) = 
h\ (Iq) .parametrized by the real parameter A G 8 (—tt/2, tt/2)) These Hamiltonians are spec- 
ified by the asymptotic s.a. boundary conditions at the origin 

V>a(p) =c[p 1/2 In {yij2p\ cos A + p 1/2 sin a] + 0(p 3/2 lnp), p -> 0, (32) 

where c is an arbitrary constant, and their domains are given by 

D hx (l ) = W- i> e ^ (Jo) (»+) , V> satisfies (E2D}, (33) 
£>£ (fa) (K+) = : V>* are a.c. in E + , z/>*, h (l ) V* G L 2 (R + )}. 

The spectrum of /ia (^o) is simple and discrete and is bounded from below, spech\ (l ) = 
{£ m , m G Z + }, where £ m are solutions of the equation u\(£ m ) = 0, 

uj x (W) = cos X[ip(a ) - 2V>(1)] - sin A, a = 1/2 - W/2 7 . (34) 

In the case of A ^ ±vr/2, the limit A — > ±7r/2 in these solutions yields the spectrum in the 
case of A = ±7r/2. 



14 



The eigenfunctions U^ m of the Hamiltonians h\ (l ) are 

U \,L = Q A,m [ui (p; 8 m ) sin A + u 3 (p; £ m ) cos A] , 
Ul (p;W) = ( 7 /2) 1/4 p 1/2 e-^ 2 / 4 $(a ,l;7P 2 /2), 

w 3 (p; WO = ^0 ln (Vt72p) + (7/2) 1/4 p 1/2 e" 7p2/4 9 M $(l/2 + p - W/2 7 , 1 + 2p; 7P 2 /2)| M=0 



Qx, 



Cu\{E v 



,1/2 



, wa(W) = sinA[^(a ) - 2^(1)] + cos A, 



27u/ A (£ m 

^ (/ ) u{% = 8 m u[% , meZ + , (35) 

they form a complete orthonormalized set in the Hilbert space £i . 

We note that spectrum and eigenfunctions in the case of A = ±7r/2 can be obtained from 
the respective formulas for the first region in the formal limit I — > 0. 

Complete spectrum and inversion formulas In the previous subsubsecs., we con- 
structed all s.a. radial Hamiltonians h e (I) associated with the s.a. differential operation 
h (I) as s.a. extensions of the symmetric operator h(l) for any / £ Z and for any any <p Q , p, 
and B. We assemble our previous results into two groups. 
For p = 0, we have 

h c (l) = h {1) (l), l^l = 0, D h(1){l) = Dl {l) (R + ), 

h(lo) = h x (k), AG§(-7r/2,7r/2) , (36) 

and the domain D hx (i Q ) is given by eq. (1331) ; 
For p > 0, we have 

h t (l) = h {1) {l), l^l a = a = 0,-l, D h(1){l) = Dl (l) (R + ) , 

h e (la) = k a (la), A a eS(-7r/2,7r/2), (37) 

and the domain Dh Xa (i a ) is given by eq. (I29I) . 

As a result, we find a family of all s.a. rotationally-invariant two-dimensional nonrela- 
tivistic Hamiltonians = M^Ti^ associated with the s.a. differential operation H 1 - (joT) 
with 5^0. Each set of possible s.a. radial Hamiltonians h e (/) generates a s.a. rotationally- 
invariant Hamiltonian in accordance with the relations (1151) and (Flo]) . As in the case of 
pure AB field where B = 0, we let E denote the spectrum points of H^~. 

It is convenient to change the indexing I, m of the spectrum points and eigenfunctions 
to I, n, as follows: 



n = nil, m) = I m ' , ~f~ . ^ _ , m G Z+, / G Z; 

' m + /, / > 

m = vain, I) = t U ' ^ ,~ _ ^ , ^ , n G Z + , / G Z, (3£ 

1 n — /, < / < 77 



and then to interchange their position, such that finally, the indices I, m are replaced by 
indices n, I. 
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When writing formulas for eigenf unctions \P n / of the operator Hj- in terms of eigen- 
functions U^ m of the operators h e (l), we has to introduce the factor l/ v / 2np e l€ ^°~ 1 ^ in 
accordance with eq. (|S]) and to make the substitution £^ m = ME n j for the corresponding 
spectrum points. 

The final result is the following. 

There is a family of s.a. two-dimensional nonrelativistic Hamiltonians parametrized 



by real parameters A*, = H 



#i = £®£- L (0©£®^i(i.)> 

1=1, 

H 1 (l) = M- 1 S l h {1) (l)Sr 1 ,l^k, 

Hl(Q = M- 1 S u h K (l a )Sl\ 

( l , n = f AG§(-7r/2,vr/2), /i = 

* \/ a ,/i>0' A * \ A a G§(-7r/2,7r/2), fi>0 ■ [6J) 

The spectrum of H^is given by (remember that E n j = M _1 8i tm i n ,i)) 

specif = {U, ez ,i^ n G Z + }} U {u, =z . {£< A *>, n G Z + }} , 

£ n ,, = 7 M" 1 [1 + 2ri + 20(Z)/x], l<n,l±U, 9(1) = j J' 

" 1 £^ /2) = 7 M~ 1 (l + 2n) 



1, I > 
/ < 



(40) 



E [ n a) = M- 1 [r a>n + 7 (a + /i)] , cu Aa (r a , n ) = 
Ei ±V2) = 7 M" 1 [1 + 2n + 26(a) /j] 



(±tt/2) . ,,-lr, .w„ , ' n 

G Z + , /i > 0, (41) 



where u\ (W) and u\ a (W) are given by the respective eqs. (|34|) and (!30|) . 

The complete set of orthonormalized eigenfunctions of consists of the functions ty n ,l(p)i I ^ 
h, and ^X(p) , 

<MP) = 7^e^^-O^W (no(p)) (42) 
where u\\l(p) are given by eqs. ( 1271) . and (we note that m(n, /*) = n) 

*iU(p) = ^e fa *»*'l7j£(p), ai = 0, 

where U\ n (p) and U\' n (p) are given by the respective eqs. fl35l) and fl3TT) . such that 

^i^(p) = EnflM, I + h, Hi(p)^ X n X = E^^SP)- 

We note, that for the case of A = ±7r/2 where / = Iq = and for the case of A a = ±7r/2 where 
I = l a = a = 0, — l,.the energy eigenvalues En and E^ a \u.d corresponding eigenfunctions 
\I/ A and ty* a are given by the respective eqs. (j4"0~|) and f|42|) extended to all values of Z. 
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The corresponding inversion formulas are 

*(p) = E E^-^Mph E $ ^**i»> 

ZeZ, Z^Z* raGZ+ 1=1*, ni 



2.2 Three-dimensional case 

In three dimensions, we start with the differential operation H (J3J). The initial symmetric 
operator H associated with H is defined in the Hilbert space Sj = L 2 {R 3 ) by 



D H = V(R 3 \R Z ), Hip = Hip, Wip G D H , 

where P(M 3 \IR Z ) is the space of smooth and compactly supported functions vanishing in a 
neighborhood of the z- axis. The domain Dh is dense in $), and the symmetricity of H is 
obvious. A s.a. quantum Hamiltonian must be defined extension of H . 

There is an evident space symmetry in the classical description of the system, the sym- 
metry with respect to rotations around the z-axis and translations along this axis, which is 
manifested as the invariance of the classical Hamiltonian under these space transformations. 
The key point in constructing a quantum description of the system is the requirement of the 
invariance of the quantum Hamiltonian under the same transformations. 

Namely, let G be the group of the above space transformations S: r — > Sr. This group 
is unitarily represented in fj: if S G G, then the corresponding operator Us is defined by 

(U s ip)(r)=iP(S- 1 r),Vipef). 

The operator H evidently commuted with Us for any S. We search only for s.a. extensions 
H e of H that also commute with Us for any S. This condition is the explicit form of the 
invariance, or symmetry, of a quantum Hamiltonian under the space transformations. As in 
classical mechanics, this symmetry allows separating the cylindrical coordinates p, tp, and 
z and reducing the three-dimensional problem to a one- dimensional radial problem. Let 
L 2 {R x R + ) denote the space of square-integrable functions with respect to the Lebesgue 
measure dp z dp onlx R + , and let V : Ylt^z L 2 (R x E + ) — > f) be the unitary operator defined 
by the relation 

2n ^p J*, 

In the two-dimensional case considered in the preceding subsection, all s.a. Hamiltoni- 
ans were represented as direct sums over I of s.a. extensions of suitable radial symmetric 



3 We remind the reader of the notion of commutativity in this case (where one of the operators, Us, is 
bounded and defined everywhere): we say that the operators H and Us commute if UsH C HUs, i.e., if 
4) G Dh, then also Usip G Dh and Us Hip — HUsip 
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operators. In the three-dimensional case, a new continuous parameter p z arises in addition 
to the discrete parameter /, and analogous representations should be written in terms of 
direct integrals rather than direct sums. More precisely, any s.a. Hamiltonian H t can be 
represented in the form 



H t = v f d Pz J2®k(i,p z )v-\ 

lei 



where h e (l,p z ) for fixed / and p z is a s.a. extension of the symmetric operator h(l,p z ) = 
h(l) +p 2 z /2m e in L 2 (IR + ), and the operator h(l) in L 2 (R + ) is given by 

D Hl) =V(R+), 

Hl)f(p) = (~d 2 p + p- 2 + P + IP 2 72) 2 - 1/4]) f(p), f G D m . (43) 

A detailed derivation of the above representation for H e will be published in the short run. 

The inversion formulas in three dimensions are obtained by the following modifications 
in the two dimensional inversion formulas: 

1) J2mez I dE — > f dp z J2iezI dE L , where E 1 - are spectrum points of two-dimensional 
s.a. Hamiltonians H^~, whereas the eigenvalues (spectrum points) E of three-dimensional 
s.a. Hamiltonians H e are E = E 1 - + p 2 z /2m ^ f dp z = dp z ^j . 

2) The contribution of discrete spectrum points of two-dimensional s.a. Hamiltonian 
have to be multiplied by f dp z . 

3) Eigenf unctions of two-dimensional s.a. Hamiltonians if^have to be multiplied by 
{2nh)~ l l 2 e iPzZ / h in order to obtain eigenfunctions of three-dimensional s.a. Hamiltonians 
H t . 

4) The extension parameters A a and A have to be replaced by functions \ a {p z ) and \{p z ). 
2.2.1 S.a. nonrelativistic Hamiltonians with AB field 

For the case of \x = 0, there is a family of s.a. three-dimensional Hamiltonians parametrized 
by a real-valued function \(p z ) (X(p z ) £ S (— 7r/2, 7r/2) , p z G M). 
The corresponding spectrum of H\( Pz ) is 

ft _ f p 2 /2m e -4M- 1 ^ e xp[2(tanA(p 2 )-C)], \X(p z )\ < vr/2 \ 

The complete set of ortho normalized generalized eigenfunctions of H\( Pz ) consists of func- 
tions *,„^(r), I ? k, and *j^(r), 

^ )E x(r) = (8n 2 h/M) ~ 1/2 e ip^+^<t>o-i)fj^ (VM£V), 

*^ Ei (r) = (8vr 2 ft (Y + 7T 2 /4) /m) jp**/*+it*<hv \j Q (Vm£>) + ^N (VME ± p) , 
A = tan X(p z ) - C - In ^v / ME x /2/t ) , 

and functions ^i^( r ), 
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jp z z/h+ie q (f>Qtp 



2M 2 



E 



M-) 



Kn 



0, \{p z ) = ±n/2 



M 



E 



■-L(-) 



p , |A(p z )|<7T/2 



^\(p z ) = -4M- 1 ^exp2(tanA(p 2 ) - C), 



y A( Pz ) 

such that 



H%^ E ±(r) = {p 2 z /2m e + E L ) ^, Pz , E x(r), £ x > 0, 

#<t!^ r ) = &/ 2m * + E± ) <t!^( r )' E± * °> 

The corresponding inversion formulas are 



/_ roc poo 

d Pz £ / ^.(S- L )*, lP „ s x(r)dE- L + / $ i „,p^ 1 )< t! £ iW^ i + $io ft <t ) ( r ) 

^.(^) = / ^W4r)*(r)dr, / ^ Z , 

<W^) = / <^ B x(r)*(r)dr, = / ^£ ) (r)^(r)rfr, 

/ |^(r)| 2 rfr = / dp, £ r ^(i^)! W + 1*^1 



For the case of p > 0, there is a family of s.a. three-dimensional Hamiltonians -£/|A a ( Pz )} 
parametrized by two real-valued functions \ a {p z ) (A a £ § (— 7r/2, vt/2) , a = 0, — 1, p z G K. 
The spectrum of H^x a ( Pz )} is 

x a = |/i + a|, A a = T(l - x a )r _1 (l + x a )tanA a (p 2 ), x a = |// + a|. 

The complete set of orthonormalized generalized eigenfunctions of H^x a ( Pz )} consists of 
functions * J)P „ B x(r), / ^ Z„, and *^ E ±(r), 

* JiPxjB ±(r) = (87r 2 n/M) _1/2 e'^/a+Mtfo-OiPj^ ( VM£» , 



q^ E± (r) = (8ir 2 hQ a ) 1/2 e *«*/ft+*«(*o-M¥> (VME^p 



(vW^o) 2 "" A a J_^ a (VME^p 



2>t n A 2 

a' 



Q a = 1 + 2 (M^/4)^ A a cos(7rx a ) + (M£/4) 2x ° A 
and functions ^M Pz) (r), 
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<:?; ) W = (2vr 2 n)- i e 



2 h\ 1 JPz z / h + i ( l a+£q<t> )<P 



M 2 E^ ( ~ ) , sin(7r^ a ) 

Xg(pz) X — K 



2-ITXa 

0, K{p z ) £ (-vr/2,0) 



M 



E ±(-) 

Aa(Pz) 



P , ^a{Pz) G (-7T/ 



■J-(-) 



A a ( Pz ) " -4M- 1 ^exp2(tanA a (p 2 ) - C), 
such that 



&%^±{t) = (p 2 z /2m e + E 1 ) % jPz>E ±(r), E L > 0, 
S <tUr) = (Pl/ 2m e + E±) <t:ix(r), ^ > C 

H<t\r) = {til*** + • 
The corresponding inversion formulas are 

//*oo 

Jez, Z^a 

a 

^ P XE L ) = J M^x(r)*(r), S 1 - > 0, Z ^ Z a , 

*J.*(^-) = / dr*^ ± (r)«r(r), & > 0, $^ = J dv^\v)^{v) 
dr\^(r)\ 2 = dp z \^ Pz (E^)\ 2 dE^ + J2^ 

J Je-rp JO „ 



, V*GL 



2/to>3\ 



2.2.2 S.a. nonrelativistic Hamiltonians with magnetic-solenoid field 

There is a family of s.a. three-dimensional Hamiltonians H\ t (p x ) parametrized by real- valued 
functions \*(p z ) (K(Pz) £ §> (— 7r/2, 7r/2), p z Gl), A* are defined by eq. ( |39l) . 
The spectrum of #a„(p») is 

spec£ A .(p,) = {p'/2m e + 2^ A *H n E Z + } U frM" 1 , 00), 

where are defined by eqs. (HOI) and (HIT) with the substitution A* — > A*(j> z ). 

The complete set of generalized orthonormalized eigenfunctions of ^A»(p s ) consists of 
functions * Pjt)ijn (r), Z 7^ Z„ and ^^(r), n G Z+, 



i)VP7> 1 r ! ■ • (44) 
where Z* are defined by eq. (139]) . m(n, Z) is given by (1551) . and uj-^p) are given by eqs. (T2j 



Vj/*(p*) C r ^ _ ip z z/h+ie(j> n <pTA3) / \ _ q 



,T r A a (p z ) / x _ 1 ip z z/h+ie(4> n -l a )tpTA2) / \ >f) 
*p„I«,nVP; - 2n^/Kp U A a ( Pz ),nlPJ' Zi > U, 
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where U\^ p ^ n {p) and U^,,{p) are given by the respective eqs. fl35|) and ([51]) with the 
substitution A* — > A*(p 2 ), such that 

^„,,„(r) = (p 2 /2m e + ^) * p „,, n (r), I ^ /*, (45) 

where 

= 7 M" 1 [1 + In + 20(Z)//], l<n, l^h, ne Z+. (46) 

We note that for A(p z ) = A a (p z ) = ±7r/2, the energy eigenvalues and corresponding eigen- 
functions \l/p zi / jn (r) are given by eqs. (I4"5l) . (1461) . and ( 1441 extended to all values of L 
The corresponding inversion formulas are 



®P*,l,n%,lMn,l)(r)+ ®P*l*^ X p lii(r) 
Z, l^U l=U 



*(r) = d Pz J2 
%^, n = J drtf p „ Jin (r)tf(r), i ^ J,, = | drtf^Jr)* (r), 

/ dr |M/(r)| 2 = [dp z J2Yl \® P ^n\\ V* G L 2 (M 3 ). 

3 S.a. Dirac Hamiltonians with magnetic-solenoid field 
3.1 Generalities 

In this section, we set c = ft = 1. Written in the form of the Schrodinger equation, the Dirac 
equation with the magnetic-solenoid field for a relativistic particle of mass m e , spin 1/2, and 
charge q = e q e ( positron or electron) .is 

r)fy ( r ) 

i^- 1 - = H*(x),x= (x°, r) , r = (x k , k = 1, 2, 3) , x° = t, 

where ^> (x) = {ip a (x), a = 1,...,4} is a four-spinor and H is the s.a. Dirac differential 
operation, the " formal Dirac Hamiltonian" , 

H = at (p — e q eA) + m e /3, 

where the vector potential A is given by fl2J), & = (7°7 fe , k = 1, 2, 3) , (3 = 7 , and 7^, \i = 
0, 1, 2, 3, are the Dirac 7 matrices. 

The space of quantum states for a particle is the Hilbert space = L 2 (M 3 ) of square- 
integrable four-spinors ^(r) with the scalar product 

(^1,^2) = J dr^f(r)^ 2 (r), dv = dx 1 dx 2 dx 3 = pdpdtpdz. 
The Hilbert space fj can be presented as 



0=1 
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Our first aim is to construct a quantum Dirac Hamiltonian that is a s.a. operator associated 
with differential operation H in this Hilbert space. A construction is based on the known 
spatial symmetry in the problem^, which allows separating the cylindric coordinates p, ip, 
and z, and on the s.a. extension theory. 

It is convenient to choose the following representation for 7 matrices: 



7 



7 



a 3 
-a 3 



7 



10 




2 > 



_°j J ) ' 7 5 = -i7°7 1 7 2 7 3 = - 



7 = 

/ 

1 



-1a 1 








ia 


S 3 = 


( 



a 3 
a 3 



Written in the cylindric coordinates, the differential operation H then becomes 



H 



Q 



cr d p + p \id v + e q (p 
ia 3 d z Q 

= e Oo + V + 7P 2 / 2 )> es0 
Q = a 1 sin <p — a 2 cos <p = — iantidiag (e i</3 , — e~ 4</3 ) 



TfipCT 3 —ia 3 d z 

o- 3 d p + p _1 (ic^ + e 9 <; 



+ p, O = [e B 0], < p < 1, 7 = e|S| > 0, 

" e* M =L 



The formal Dirac Hamiltonian commutes with the s.a. differential operations 



Pz 



id z , S z = 7 5 (7 



" ///, ; p :: ) 



Jz = L z + -T 3 = -id v + -T 3 = diag (j„ j z ) , j z = -id v + cr 3 /2. 

We pass to the p z representation for four-spinors, \I/(r) — > ty(p z , p), 
1 , , , = , , 1 



In this representation, the operation J 2 is the same, while the differential operation H and 
the operation S z respectively become 



H^H(p z , 

S z -> £ 2 (Pa 







cr 3 (9 p + p 1 ( ic? v + 



a 3 p z 



m e L p z 



+ m e a 3 a 3 p z 

Q 



a 3 d p + p 1 (id v + e 9V 



— m P a" 



m e 1 p z 



We decompose the four-spinor ty(p z ,p) for a fixed p z into two orthogonal components 
that are the eigenvectors for the spin matrix S z (p z ): 



By the spatial symmetry, we mean the invariance under rotations around the solenoid axis and under 
the translations along this axis. 
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The space of four-spinors ^f(p z ,p) with fixed p z is the direct orthogonal sum of two 
eigenspaces of S z (p z ), 

#(p z , p) = *i(p z , p) + *_i(p z , p), 



where 



M ~ 



S z (Pz) ^ s (Pz, p) = s—^ s ( Pz , p), M = y/m 2 e +p>, s = ±1, 

77l e 

*.<*. p) = (^^) ( ft K + jv/) -1 Xl ) = P) ® «i(P.), 

*-.<*.,*> = (^) V2 ( -"«( m - + f)" 1 ^ ) =X-,fc,p)»e_ lW , 



e s(Pz)? s = ±1, are two orthonormalized two-spinors, 

(2M)- 1 / 2 (m e + M) 1/2 
p z [2M(m e + M)] 



= ( [ 'ZL ) , *-i(P.) = -^i(P,), (47) 



and Xs(Pz) P) are some doublets. We thus obtain a one-to-one correspondence between four- 
spinors ^(r) and pairs of doublets x s (Pz, p), 

tf(r) * s (p„ p) Xs fe, P), 

such that 1 1^| | 2 = EJIXsll 2 = Yls I d Pzdpxf(Pz, p)x s (Pz, p)- 

The differential operations ij and J 2 induce the differential operations h and j z in the 
space of doublets x s (Pz, p) : 

H (p z ) ^ s = h (s, p z ) x s ® e s, Jz (Pz) = JzXs ® e s , 



h(s,p z )=Q a 3 d p + p 1 (id ip + e q ^j 



sMa c 



The s.a. operator j z associated with the differential operation j z has a discrete spectrum, 
its eigenvalues are all half-integers labelled here by integers I as e(0 o — Z+l/2). It is convenient 
to represent vectors Ci(f) = £i(Pz,P,<p) of the corresponding eigenspaces, 

j^) = [e(0o-* + V2)]6(¥>), / GZ, 



as 



= (2 7 r)- 1 / 2 e l ^o-^i/ 2 )- 3 / 2 ]^ = Sl (ip)-±=F(l,p z , p), 

Sfa) = e^°- m/2 ^antidiag (ie 1 ^ 2 , -e^ 2 ) , S+(<^) Sfa) = I, (48) 

where $i = $i(p z ,p) and F(l,p z ,p) are arbitrary doublets independent of ip. 

The space of each of doublets x s (Pz, p) is a direct orthogonal sum of the eigenspaces of 
the operator j z , which means that the doublets allows the representations 

X s (Pz,p) = ^2-^=Si((p)F(s,l,p z ,p), 
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the factor 1 / y / 2irp is introduced for further convenience. 

The operation h (s,p z ) induces an operation h (s, I) ("radial Hamiltonian" depending on 
parameter p z as well) in the space of doublets F, 

h (s, I) = ia 2 d p + e(7p/2 + p^xfia 1 - sMa 3 , (49) 

where x\ = I + p, — 1/2. 

In the Hilbert space JL 2 (R + ) = L 2 (R + ) © L 2 (R + ) of doublets F(p) (with p 2 fixed), we 
define the initial symmetric radial Hamiltonian h (s, I) associated with the s.a. differential 
operation h (s, I) by 

h{Sj) -\h(s,l)F(p) = h(s,l)F(p). ^ 

3.2 Solutions of radial equations 

I. We first consider the homogeneous equation 

[h(s,l)-W]F{p) = (51) 

and some of its useful solutions. 

We let / and g denote the respective upper and lower components of doublets F, 
F = {f/g) ■ Then eq.f lSTj) is equivalent to the system of radial equations for the doublet 
components 

/' - e( 7 p/2 + p- 1 ^)/ + (W- sM)g = 0, 

g' + e( 7 p/2 + p^x^g - (W + sM)f = 0, (52) 

where the prime denotes the derivatives with respect to p. 
We start with the case of e = 1. 

The system ( l52p can be reduced to second-order differential equations both for / and g. 
For example, we have the following system equivalent to (1521) : 



(7P/2) H - w + 7 X/ + 



/ = 0, 



p 2 ~ ' ' V'" ' 2 / 

g=(W- sMY 1 [-/' + ( 7 p/2 + p- 1 ^)/] , «; = - M 2 . (53) 

By the substitution 

f{p) = z a / 2 e- z ' 2 p{z), z = 7 p 2 /2, a = 1/2 ± (x, - 1/2), 

we reduce the first equation (|53|) to the equation for p(z) that is the equation for confluent 
hypergeometric functions, 

zd 2 p + (/3 - z)0*p - ap = 0, /3 = a + 1/2, (54) 
a = a/2 + Xi/2 + 1/2 - w/2 7 . 
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Known solutions of eq. (15"4"j) allows obtaining solutions of eqs. ( l5"Tj) . 

In what follows, we use the following solutions F\(p\ s, W), F 2 (p; s, W), and F 3 (p; s, W) 
of equation (lo"Tj) : 



l/2-Z- M -z/2 



- (2/?!)- 1 (W - sM)p$(ai + 1, /?! + 1; z) 



F 2 = p^-^e 



(2/? 2 )" 1 (W + sM)p$(a 2 , & + !;;*) 



F 3 -p< e ( ^ (ai;/3i . z) ) - ^1 - ^F 2 , (55) 



where 

(3 X — 1 — I — p, «i = — w/2j, (3 2 = I + p, 0*2 = I + p — w/2j, 

All the solutions Fx, F 2 , and F3 are real-entire in VF. 

The solutions (1551) have the following asymptotic behavior at the origin and at infinity. 
As p — > 0, we have: 

F 1 = pV--^-(2/5 i r 1 (^-,M)p\ (l + 0(p2))) 
F -^ 1/2 ((2^(^ + S M)p)( 1 + °^ 

/, - nr -;- u)rMi) x J a + °oP-4 frio, 

2(7/2) ^(oti + l) (1 + 0(p 2 In p)), Z = 0, p = 

„_W 1/2 ^ x f (l + 0(p 2 )), />! 

^ 3 ~i> 2 ) p x \ (i + o(p^)),/ = o, p>o (56) 

where F 3 = {fo/gz) . 
As p — > 00, we have: 

F 2 = (7/2 ^^ 2) p" +2Ql e-/ 2 (1/ (7P)- 1 (W + sM)) (1 + 0(p- 2 )), 
F 3 = ( 7 /2)" ai f-^e-'/ 2 ((7P)" 1 (W - sM) /l) (1 + 0(p~ 2 )). 

We define the Wronskian Wr(F, F) of two doublets F = (/ ' / g) and F = ^f/gj by 

Wr(F,F) =j~g-g] = iFo 2 F. 
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If (h - W)F ={h- W)F = 0, then Wr(F, F) = C = const. Solutions F and F are linearly 
independent iff C ^ 0. It is easy to see that Wr(F 1; F 2 ) = — 1. 

If Im W > 0, the solutions Fx, F 2 , and F 4 are pairwise linearly independent, 

Wr(Fi,F 3 ) =ux{W), Wr(F 2 ,F 3 ) =u 2 (W). 

Taking the asymptotics of the linearly independent solutions Fx and F 3 into account, we 
obtain that there are no square integrable solutions of eq. flSlH with ImW ^ and |/| > 1 
or I = 0, p = 0. This implies that in these cases, the deficiency indices of h (s, I) are zero. 
In the case of I = 0, \i > 0, the solution F 3 is square integrable, which implies that the 
deficiency indices of h (s, 0) are equal to (1,1). 

For any I and [i, the asymptotic behavior of any solution F of eq. (I5TT) at the origin, as 
p — > 0, is not more singular than p - '^', -F(p) = 0{p^^). 

II. We now consider the inhomogeneous equation 

(h (s, I) - W)F(p) = V(p), W G L 2 (M + ). 
Its general solution allows the representations 
F(p)=CxF d (p; W) + c 2 F 3 (p; W) 



F d (p;W) / F 3 (r;iy)^(r)dr + F 3 (p;W0 / F d {r; W)V{r)dr 



u d = Wr(F d , F 3 ), d=l,2, I < for d = 1, I > 1 for d = 2. 



(57) 



A simple estimate of the integral terms in the r.h.s. of (joTj) using the Cauchy-Bunyakovskii 
inequality shows that they are bounded as p —>■ oo. It follows that F 6 L 2 (R + ) implies ci = 0. 

For \xi\> 1/2, an evaluation shows that as p — > 0, the integral terms are of the order of 
0(p 1 ^ 2 ) (up to the factor lnp for \>c\\ = 1/2) . In this case, F G L 2 (IR + ) implies c 2 = 0, and 
we find 



F(p) = u 



-i 



F d {p;W) / F 3 (r;W)q(r)dr + F 3 (p;W) / F d (r; W)^(r)dr 



(58) 



For < 1/2, the doublet F 3 (p; W) is square-integrable, and a solution F(p) G L 2 (R + ) 
allows the representation 



F(p) = bu^Fxip; W) + c 2 F 3 (p; W) 



^3(P;^) ri r i(r;W / )^(r)c/r-F 1 (p;I^) F 3 (r;W)ty[ 
Jo Jo 



r)dr 



F(p) = bio^ L Fx(p; W) + c 2 F 3 (p; W) + 0(p 1/2 ), p -> 



(59) 



where 



F 3 (r;W)m(r)dr . 



We use representations. fl57|) - fl59|) to determine the Green functions for s.a. radial Hamil- 
tonians. 
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3.3 S.a. radial Hamiltonians 
3.3.1 Generalities 

We proceed to constructing s.a. radial Hamiltonians h e (s,l) in the Hilbert space L 2 (M + ) of 
doublets as s.a. extensions of the initial symmetric radial operators h (s, I) fl50l) associated 
with the differential operations h (s, I) (|49[) and analyze the corresponding spectral problems. 

The action of all of the following operators associated with the differential operations 
h (s, I) is given by h (s, I), therefore we cite only their domains. 

We begin with the adjoint h + (s, I) of the initial symmetric operator h (s, I). Its domain 
Dh+ is the natural domain for h (s, I) , 

D h+ = D h(s,l) = { F * (P) : F * aX - 111 M +' F *> h ( S > F * e L2 ( M +)} • 

The quadratic asymmetry form (F#) of h + (s, I) is expressed in terms of the local 
quadratic form 

[F*, F.] (p) = W)J{P) - W)9(P), F * = (f/g) , 

as follows: 

A h+ (F*) = (f^+F*) - (h + F*,F^ = - [F m ,F m ] . 

We can prove that linip^oo F* (p) = for any F* G D~(R + ). Indeed, because F* and 
h (s, I) F* are square integrable at infinity, the combination 

Fi - (tp/2) a 3 F* = -ia 2 [h (s, Z) F* - (x l /p)a 1 F, + sMa 3 F # ] 

is also square-integrable at infinity. This implies that / and /' — (jp/2) f, together with g 
and g' + (jp/2)g, are square-integrable at infinity. We consider the identity 



\f(p)\ = / [df(r)f(r) + f(r)df(r)]dr + 1 r\f(r)\ 2 dr + \f(a)\\ d = d p - jp/2. 

J a J a 

The r.h.s. of this identity has a limit (finite or infinite) as p —>■ oo. Therefore, \ f(p) \ also has 
a limit as p — > oo. This limit has to be zero because /(p) is square-integrable at infinity. In 
the same way, we can verify that g(p) — > as p — ► oo. 

To analyze the behavior of F* at the origin, we consider the relation 

V = h(s, l) F m , #, F* G L 2 (R + ), (60) 

or 

/' - (TP/ 2 + P -1 ^) / = -X2, 9 + (TP/ 2 + P" 1 ^) £ = Xi, 
X = (X1/X2) = ^ + sMa^F* G L 2 (M + ), 

as an equation for F* at a given \. The general solution of these equations allows the 
representation 



f{p) = P^ p2/4 



ci+ / r- >Cl e^ r2/4 x 2 (r)dr 



g{p) = p 



-^io-7P 2 /4 



c 2 + [ P r^e^ 2/4 Xl (r)dr 



(61) 
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It turns out that the asymptotic behavior of the functions / and g at the origin crucially 
depends on the value of I. Therefore, our exposition is naturally divided into subsections 
related to the corresponding regions. We distinguish three regions of 



3.3.2 First region: k x < -1/2 
In this region, we have 

/ < 



-1, p > 
0, p = ' 

The representation ( |6T|) allows estimating an asymptotic behavior of doublets G 



D^ sl j at the origin for the first region. 



f(p) = p-Merf* 



p 



~c x - \ r^e^ /4 x 2 (r)dr 



o 



c^I+OCp 1 / 2 ), p^O, 



oo 

,2, 



ci = a+ / r^e-^ /4 x 2 (r)dr. 



o 



The condition / G L 2 implies c\ = 0, and therefore, f(p) = 0(p 1 ^ 2 ) as p — ► 0. As to 
g(p), we find 

/ oip 1 / 2 ), ^ < -i/2 

9{P) - \ 0{p l / 2 \np), = -1/2 (/ = 0,/i = 0) ' P ^ U - 

We thus obtain that (p) — > as p — > 0, which implies that A fe + (F*) = 0, VF* G 
.D?^ (R+) • This means that the deficiency indices of each of the symmetric operators h (s, I) 

in the first region are zero. Therefore, there exists only one s.a. extension h t (s,l,p z ) = 
(s,l) = h + (s,l) of h(s,l), i.e., a unique s.a. radial Hamiltonian with given s and I, its 
domain is the natural domain, i ^ ^ = Di, „ (R+) • 

The representation (158|) with d = 1 implies that the Green function for the s.a. Hamil- 
tonian h(i) (s, I) is given by 

C(o o>- W) - -J— l W) ® Fl{p ' ] W ^P>P' (62) 

G ^ P ' W) -Z^W)\ p<p' • (62) 



Unfortunately, we can not use representation (I55j) for F3 as a sum of two terms directly for 
all values of p because the both are singular at p = (although the sum is not). 

To cover the total range of p, we use another representation for F3.. 

We let F d i(p; W) denote the functions F d (p; W), d = 1, 2, 3, with a fixed I and represent 
F31 as 

F 3l = ui[AuFu + F 4l ], An = A U (W) = Q 1 (W) - r(f3 2 )P u (W), 

uj 2 {W) 



Fy = Fu(p; W) = T(j3 2 )P il {W)F Il ip; W) - F 2l (p; W), Sl^W) 
Pu{W) 



ui(wy 

(W + sM)( 1 /2fT(a 1 ) 



2|/|!r(a 1 -|/|) 
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Using the relation (see |24|) 

1 x n+1 T(a +n + l) 

lim ^$(«, /3; x) = — J $(a + n + l,n + 2;x), (63) 

we can verify that 

r-\[3 2 )F 2l (p;W)\^ = P ll (W)F ll (p;W)\ tl=0 . 

Taking the latter relation into account, it is easy to see that in the first region, An and F& 
are finite for p > 0, as well as u)\ and Fu, and also that Pu(E) and F^^p; E) are real. 
The Green function is then represented as 

G(p, p'; W) = A ll (W)F ll (p; W) ® F u (p'; W) 

+ \ F ll (p;W)®F 4l (p>;W), p < pf ' l ° J 

for all /i > 0. 

We choose the guiding functional $i(-F; IU) for the s.a. operator hm (s, I) in the form 

Fi(p; W)F(p), F(p) GD = D r (M + ) n 

It is easy to prove that the guiding functional is simple. It follows that the spectrum of 
h(X) (s, I) is simple. 

Using representation (1M|) for the Green function , we obtain that the derivative cr'(E) = 
[irFi (p; W)} -1 ImG(p, p; E + iO) of the spectral function is given by 

a'(E) = n- 1 lmA ll (E + iO). (65) 

It is easy to prove that Im An{E + iO) is continuous in p for p > 0, such that it is sufficient 
to find cr'(E) only for the case of p > where eq. (1651) is more simple, 



(W + a M)( 7 /2)-*r(/3 a 
0" (£J = — 



2 7 rr(/3 1 )r( 0;2 ) 



Imr(ai)| w=B+i0 . (66) 



W=.E+iO • 

W=E 



It is easy to see that cr'(E) may differ from zero only at the points Ek defined by the 
relation a\ = —k (T(ai) = oo), or M 2 — E\ = —2^k, which yields 

E k = ±M k , M fe = ^/M 2 + M = M, k e Z + . 

The presence of the factor (E + sM) in the r.h.s. of (I66I) implies that the points E = —sM = 
—sM do not belong to the spectrum of h(i) (s, I) . In what follows it is convenient to change 
the numeration of the spectrum points. Introduce an index n(s): 

n(s) G Z{s) = K(s)} , a = ±, n+(s) G | ^ S = ~_\ , n_(s) G j ^'/j^ • 
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Then we have 

E k = ±M k =► E <a) = (rM| n(s) |, n(s) G Z(s). 

Finally, we obtain 

a'{E)= J2 Ql(s)S(E-E n(s) ), 

n(s)e2(s) 



; ( 7 /2)^r(/j 1 + |n( g )|)(l + g ME fc - 1 ) 
= \/ In^llPCft) ' /?1 = 1 + 1/1 " 

Thus, the spectrum of the s.a. Hamiltonian /i^) (s, Z) is simple and discrete, specfyi) (s, I) = 
{E n(s) , n(s)EZ(s)}. 
The eigenvectors 

i i 

U n (s) = U n{s) (s,l,p z ;p) = Q n (s)F 1 ( y p]E n{s) ), n(s) G (67) 

of the Hamiltonian /in) (s, I) form a complete orthonormalized set in the space L 2 (M + ) of 
doublets F (p). 

3.3.3 Second region: Hi > 1/2 

In this region, we have I > 1. 

The representation (l6Tj) yields the following estimates for an asymptotic behavior of 
doublets G „ at the origin for the second region: 



f(p) 



0{pV% x, > 1/2 

©(p 1 / 2 In p), x ; = 1/2 , p->0. 



= ^(p 1 / 2 ) 

It follows that F* (p) -> as p -> 0, which implies that A h + (F*) = 0, VF* G ;) (JR+) . 

This means that the deficiency indices of each of the symmetric operators h (s, /) in the 
second region are also zero. Therefore, there exists only one s.a. extension h e (s,l,p z ) = 
h(2) (s,l) = h + (s,l) of h (s,l), i.e., a unique s.a. radial Hamiltonian with given s and I, its 
domain is the natural domain, D~ h 2 ( s /) = Dh s n (M+) ■ 

The representation (1351) with d = 2 implies that the Green function for the s.a. Hamil- 
tonian hp) (s, I) is given by 



G(p,p';W)=co 2 l (W) 



F n {p-W)®F 2l {p>-W), p> p' 
F 2l {p-W)®F n {p'-W), p<p' ■ 

Again, the representation (155]) for F3 as a sum of two terms is not applicable directly for 
p = We therefore use the following representation for F3: 



F 3l = co 2 [F 5l - A 2l F 2l ], A 2l = A 2l {W) = tt 2 {W) + T{(3 x )P 2l {W), 



F 5l (p; W) = F u (p; W) + T((3 1 )P 2l (W)F 2l (p; W), Q 2 (W) 



u 2 (wy 



p (w) = (^- a M)( 7 /2) t - 1 r(q 1 + 

2 2(z-i)ir( ai + i) 
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Using relation (1631) . we can verify that 

r-^jFufaW)^ = - P 2l {W)F 2l {p-W\ =Q . 

Taking the latter relation into account, it is easy to see that A 2 i and F 5 i are finite for p > 0, 
as well as u 2 and F 2 i, in the second region, and P 2 i(E) and F 5 i(p; E) are real. 
The Green function is then represented as 

G(p, p'; W) = -A 2l (W)F 2l (p; W) ® F 2l (p'; W) 

f F 5l (p;W)®F 2l (p>;W), p > p' 
+ 1 F 2l (p;W)®F 5l (p';W), p < p' ' [ ™> 

for all p > 0. 

We choose the guiding functional Q 2 (F; W) for the s.a. operator h( 2 \ (s, I) in the form 

POO 

$ 2 (F; W) = / F 2 (p; W)F(p), F(p) E D = D r (R+) n D 



h(2)(s.')' 

It is easy to prove that the guiding functional is simple. It follows that the spectrum of 
h(2) (s, I) is simple. 

Using representation (168]) for the Green function, we obtain that the derivative cr'(E) of 
the spectral function is given by 

cr'(E) = -n- 1 Im A 2l (E + zO). (69) 

It is easy to prove that lmA 2 i(E + iO) is continuous in p for p > 0, such that it is sufficient 
to find cr'(E) only for the case of p > where eq. ( 1691) is more simple, 



a (£) = 



7r 7 r(/3 2 )r(i + ttl 



ImT (a 2 )|. 



W=E 



It is easy to see that cr'(E) may differ from zero only at the points E k defined by the 
relation a 2 = —k (T (a 2 ) = oo) or 

M 2 -E 2 k + 2j(l + p) = -2-ik , k e Z + , 

which yields 

E k = ±^M 2 + 2 1 (k + l + p) = ±M k+l+ ^ k e Z+ . 

All the points E k are the spectrum points. 

It is convenient to change indexing k for n(s), 

E k £ n(s) = aM| n(s) | +M , {n(s) G |n(s)| > /} (n CT (s) = cr(A; + 1), k G Z+ ) . 

We finally obtain that 

a'(E)= J2 Ql(sAE-E n ), 
nEZ,\n\>l 



Qn(s) 



f ( 7 /2)^r(|n(.)| +/,)(! - g ME-^; 
(|n( S )|-0!r 2 (/ + /i) 
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So, the spectrum of the s.a. Hamiltonian Km) (s, I) is simple and discrete, spech( 2 ) {s, I) = 
{E n{s) , n(s)eZ, \n(s)\ >l}. 
The eigenvectors 

ii ii 

U n (s) = U n{s) (s,l,p z ;p) = Qn(s)F 2 (p;E n{s) ), n(s) E Z(s), (70) 

of the Hamiltonian h( 2 ) {s, I) form a complete orthonormalized set in the space L 2 (M + ) of 
doublets F (p). 

3.3.4 Third region: |x z | < 1/2 

In this region, we have Z = Z = 0, and x; reduces to x = p — 1/2, p > 0. 

The repre 
at the origin: 



The representation ( 1611) yields the following asymptotic behavior of doublets F* E Dh ; , 



f(p) = c 1 (m e p)^ +0{p i/2 ) 
It follows that 

A h+ (F m ) = ~ 2 c x - clc 2 , m E D* h(sM (R+) . 

Up to the factor i, the r.h.s. is a quadratic form in c\ and c 2 with the inertia indices (1, 1), 
which implies that the deficiency indices of the initial symmetric operator h (s, Iq) are (1, 1). 
The additional asymptotic boundary conditions 

with a fixed A E S (— tt/2, it/2) (note that A depend on s and p z , A = X(s,p z )) on doublets 
define a maximum subspace in ^ (M+) where = 0. This subspace is the domain of 

a s.a. operator that is a s.a. extension of h (s, Iq). 

We thus obtain that there exists a one-parameter U(l) family of s.a. radial Hamiltoni- 
ans h t (s,lo,p z ) = h\(s,lo) parametrized by the real parameter A E S (— tt/2, tt/2). These 
Hamiltonians are specified by asymptotic s.a. boundary conditions (ITT]) , and their domains 
are given by 

Dk xi s,i ) = {F(p) : F(p) E Dl x(S)lo) (R + ) , F satisfies dZTj} . (72) 



According to representation ( 1591) . which certainly holds for the doublets F belonging to 
Dh x (s,i ), and (1561) . the asymptotic behavior of F at the origin is given by 



On the other hand, F satisfies boundary conditions (ITTj) . whence it follows that there must 
be 

b cos A _ 2 
c 2 = , c<j(A) = uj 2 cos A + m e ^cui sin A. (73) 

^l^(A) 
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Then representation (1591) for F with c 2 given by (1731) implies that the Green function for 
the s.a. Hamiltonian h x (s, Iq) is given by 

G(p, p'- W) = Q- 1 (W)F [x) (p; W) ® F {x) (p'; W) 
f F {X) (p-W)®F {x) (p'-Wl p>p> 

\ F {x) (p;W)®F {x) (p';W), p < p' ' 1 1 

where 

F {x) (p; W) = m-^Fxip; W) sin A + m^F 2 (p; W) cos A, 
F (A) (p; W) = m'^F^p- W) cos A - m^F 2 (p; W) sin A, 

W = ^7^> <°^3 = U>(A)F(A) + W (A) F (A)| 

uj {x) (W) 

u)(A) (W) = co> 2 sin A — m~ 2>ca uji cos A. 

We note that the doublets F^(p; W) and F( X ){p; W) are real-entire in W, and the doublet 
F( X ){p;W) satisfies asymptotic s.a. boundary conditions (J7T|) . 

We choose the guiding functional &\{F; W) for the s.a. operator h x (s, Iq) in the form 

P CO 

$ {A) (F; W) = / F (A) (p; W)F(p), F(p) G © = D r (M + ) n £> &aWo) . 
Jo 

It is easy to prove that the guiding functional is simple. It follows that the spectrum of 
h x (s, I) is simple. 

Using the representation (1731) for the Green function, we obtain that the derivative cr'(E) 
of the spectral function is given by 

a'{E) = ir-Hmtt-^E + iO). 

Because Q(E) is real, cr'(E) differs from zero only at the zero points E k of the function 
tt(E), tt(E k ) = 0, and we find. 

a'(E) = J2Ql S ( E - E k\ Qk = l-V(E)]- 1/2 , 
k 

spech x (s, l ) = {E k } , k e Z. 

The eigenvectors 

111 111 

U k = U k (\,s,p z ;p) = Q k F {x) (p;E k ), k E Z, (75) 

of the Hamiltonian h x (s,Iq) form a complete and orthonormalized set in the space L 2 (M + ) 
of doublets F (p) . 

For A = and A = ± 7r/2, we can evaluate the spectrum explicitly. 

a) Let A = 0, and we consider the s.a. Hamiltonian ho (s, Iq). We have 

F (0) (p; W) = m?F 2 {p- W), Q(W) = - e ^ ( ^ } > 
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and find 

a'(F) = - 2m ^(7/2f 2 r(ft) 
1 ' n(W + sM)r(ai)r(/3 2 ) 



ImT (a 2 )\ w=E+i0 . 



W=E 

As in the second region, a'(E) differs from zero only at the points E k , defined by the relation 
«2 = —k (T (a 2 ) = oo) or 

M 2 - El + 2 7 + p = -2jk , k G Z+ , 

which yields 

E k = ±M k+ ^, k e Z+ . 

All the points E k are the spectrum points. 
It is convenient to introduce an index n, 

n G Z = {n a G <tZ + , a — ±} , 

where n + = and n_ = are considered different elements of the set Z. Then we can write 

E k ^E n = aM H+lx , neZ. (76) 

.We finally obtain that 



a'(E) = Y,m; 2 ^Q 2 n 5(E-E n ), Q n 



( 7 /2)" r( 


n 


+ /i)(l- S MEn 1 ) 




n|!P(^) 



neZ 

So, the spectrum of the s.a. Hamiltonians ho (s, /o) is simple and discrete, spec/io (s, /o) = 

m in 

{E n , n E Z} . The eigenvectors U n = U n (0, s, l,p z ; p) = Q n F 2 (p; E n ), n E Z of the Hamilto- 
nian h (s, l ) form a complete orthonormalized set in the space L 2 (IR + ) of doublets F (p). 

We note that the spectrum, spectral function and eigenfunctions of ho (s, Iq) are obtained 
from the respective expressions for the second region, x; > 1/2, by the substitution / = 0. 
We also note that for p > 1/2, the function F( )(p; W) = m£°F 2 (p; W) is minimally singular 
at the origin among the functions F^(p; W), in fact, it is nonsingular. 

b) Let A = 7r/2, which is equivalent to A = —n/2, and we consider the s.a. Hamiltonian 
h-n/2 (s,lo). We have 

F (7r/2) (p; W) = m^F^p- W), Q(W) = } 

and find 

™ 2 *oV(R»\(W 4- s/l/fi 

ImT {oi\)\ w=E+i() . (77) 



m 2 e ^T((3 2 )(W + sM) 

( E ) = i ; ^ 



27r( 7 /2)^r(/3 1 )r(a 2 ) 

As in the first region, a'(E) differs from zero only at the points E k) defined by the relation 
cti = —k (r(ai) = oo), or 

j[/r2 772 

£ = — fc, E fc = ±a/M 2 + 27 A;, fc G Z+ . 

27 
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The presence of the factor (E + sM) in the r.h.s. of (1771) implies that the points E = —sM = 
—sM do not belong to the spectrum of h w / 2 (s, lo)- We change the indexing of the spectrum 
points : 

E k =>E n = <TM M ,neZ(s). (78) 
E k = (signAOMjjfci, \k\ > 1 E = sM ; k E Z . 
We finally obtain that 



a'(E) = J2^°Ql5(E-E n ), Q n 



r(|n| 




/i) (1 + sMSn 1 ) 


( 7 /2)* 


n 


!r 2 (i-/i) 



fcez 

So, the spectrum of the s.a. Hamiltonian (s, lo) is simple and discrete, spec/z 7r / 2 (s, l ) = 
{E n , neZ(s)}. 

in in 

The eigenvectors U n = U n (tt/2, s, lo,p z ', p) — QnFi(p; E n ), n G Z(s), of the Hamiltonian 
h-K/2 (s, lo) form a complete orthonormalized set in the space L 2 (R + ) of doublets F (p). 

We note that the spectrum, spectral function, and eigenf unctions of h n /2 (s, lo) can be 
obtained from the respective expressions for the first region, >c\ < —1/2, by the substitution 
I = 0. We also note that for \i < 1/2, the function F^/ 2 )(p] W) = m~* <0 F\{p\ W) is minimally 
singular at the origin among the functions -F(a)(p; W); in fact, it is nonsingular. 



3.4 Complete spectrum and inversion formulas for Dirac spinors 

In the previous subsubsecs., we constructed all s.a. radial Hamiltonians h e (s,l,p z ) as s.a. 
extensions of the symmetric operators h(s,l,p z ) for any s, I, and p z and for any values of 
O , /i, and 7. The total s.a. operators H e associated with the Dirac differential operation H 
in the Hilbert space f) = L 2 (IR 3 ) of Dirac spinors are constructed from the sets of h t (s, l,p z ) 
by means of a procedure of "direct summation over s, I and direct integration over p z \ 
Each set of possible s.a. radial Hamiltonians h e (s, l,p z ) generates a spatially invariant^ s.a. 
Hamiltonian H e . Namely, let G be the group of the above space transformations S: r — > Sr. 
This group is unitarily represented in fj: if S G G, then the corresponding operator Us is 
defined by 

(U s j>) (r) = e- ies3 /V(S- 1 r),V^G^, 

where 9 is the rotation angle of vector p around the z-axis. The operator H evidently 
commutes with Us for any S. We search only for s.a. extensions H z of H that also commutes 
with Us for any S. This condition is the explicit form of the invariance, or symmetry, of 
a quantum Hamiltonian under the space transformations As in classical mechanics, this 
symmetry allows separating the cylindrical coordinates p, tp, and z and reducing the three- 
dimensional problem to a one-dimensional radial problem. Let V is the unitary operator 
defined by the relation 

(yf)(j,,<p,z) = -L- I dp z Y,^ z [Si^)F{s,l,p z ,p)]®e s {p z ), 

V P -'Rz law. 



5 I.e., invariant under rotations aroud the z axis and under translations along the z axis. 



35 



where Si(<p) and e s {p z ) are given by the respective (j!71) and (flgj) . Similarly to the consider- 
ations in subsec. 2 and 3 of sec. 2, it is natural to expect that any s.a. Hamiltonian H t can 
be represented in the form of the type 



H*=V j d Pz J2J2h e (s,l,p z )V~\ 

JRz s= ±i i €Z 



where h e (s,l,p z ) for fixed s, I, and p z is s.a. extension of symmetric operator h(s,l,p z ), 

\ h(s,l,p z ) = h(s,l,p z )F(s,l,p z ,p), VF(s,l,p z ,p) G D h(s ^ Pz) , 

acting in the Hilbert space L 2 (IR + , dp) of the functions f(p,l,p z ) with the scalar product 
(F(s,l,p z ),G(s,l,p z )) = f R+ F(s,l,p z ,p)G(s,l,p z ,p)dp, h(s,l,p z ) is given by eq. (|49]). A 

correct expression for H e is 



J R z s=±l igZ 



A detailed exposition of the procedure will be published in the short run. 

The inversion formulas in the Hilbert space S) are correspondingly obtained from the 
known radial inversion formulas by a procedure of summation over s, I and integration over 
p z . And we now must consider the extension parameter A a function of s and p z , A = X(s,p z ). 
In what follows, j dp z means dp z . 

For p = 0, there is a unique s.a. Dirac Hamiltonian H e . 
The spectrum of H t is 

specif = (— oo, — m e ] U [m e , oo). 

A complete set of generalized eigenfunctions of the s.a. Hamiltonian H t is the set 

{* S)P .,«( s y(r), n(s) G Z{s), I < |n(s)|}, 

*.j>,rtfi)A T ) = l^^ iPzZ Si{^>)F n{s ){s,l,p z ] p) ® e s (p z ), 



F n{s) (s,l,p z ;p) 



U n (s) (s,l,p z ;p) , I < 0, 



n 



U n (s) (s,l,p z ;p) , 1 < I < \n(s) 



i ii 

the doublets C/ n ( s ) (s,l,p z ;p) and C/ n ( s ) (s,l,p z ;p) are given by the respective (1671) and ([70 
such that 

H^ s ,p z ,n(s),l( r ) — -^s,p z ,n( S ),«^s,p z ,n(s),/(r), 

where 

E s , Pz ,n(s),i = ay/ml+ p2 + 2 7 |n(s)|, n(s) G Z(s), Z < |n(s)|. 
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The corresponding inversion formulas are 

(r) = dp z J2 E E $ S , Pzi n( s) ,/^,, Pz ,n( s ),/(r), 
J s =±ln(s)eZ(s) l<\n(s)\ 

®s,p z As),i = J * S) p„n( s ),Kr)*(r)dr, 
I |*(r)| 2 dr = [d Pz J2 E E l*-^..»l 2 ' W G L ' ( R3 ) • 

J J s=±ln(s)eZ(s)l<\n(s)\ 

For > 0, there is a family of s.a. Dirac Hamiltonians H^\r SjPz ^ parametrized by two 
real-valued functions X(s,p z ), A G 8 (— n/2, ir/2) , s = ±1. 
The spectrum of H{\( s , Pz )} is 

specf/"{ A ( s ,p z )} = (-oo, -m e ] U [m e , oo). 

A complete set of generalized eigenfunctions of the s.a. Hamiltonian Hr\r S}P \\ is the set 
n(s) G Z(s), I < \n(s)\, I ? 0} U {^^(r), k G z}, 

^s,p z ,n( s ),K r ) = ^-^e 4p2Z 5' i (v3)F n(s) (s,/,p z ;p) ® e s (p 2 ), 



F n (s,l,p z ;p) 



i 

U n (s,l,p z ;p) , / < -1, 
ii 

U n (s,l,p g ;p) , 1 < I < |n(s) 



\c ^ 1 m 

^i ( ;Si( r ) = ^p eWzZS ^)U k (X(s,p z ), s,p z ; p) ® e s (p z ), 

in 

where U n (X(s,p z ), s,p z ; p) are given by (!75l) with the substitution A — > X(s,p z ), such that 



where, 



9(1) 



and 
where 



Es, PzA s U = a^ml+pl + 2 1 [\n{s)\+6{l)l I < |n(s)|, I ± 0, 

0, I < 

1, Z > 1 ' 



pA(^) . /x f a n n Q , A w) _ cosX + a(W) sin A 

= 2m^- 1 ( 7 /2) 1 ^r( / i)r(l- / i-^/2 7 ) 
1 J (W + sM) r (1 - n) r (-m/2 7 ) 
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We recall that for X(s,p z ) = and X(s,p z ) 
found explicitly, see the respective eqs. ( 1751) and 
The corresponding inversion formulas are 



±7r/2, the eigenvalues E^, S Z Z \ can be 



(r) = d Pz J2 

s=±l 

J ^^M^(r)rfr, / ^ 0, $ s , Pz , Mo = | <^,l (rMr)dr 



^ ^ $s, Pz ,n( S )^s, Pz ,n( S ),z(r) + $ 

n(s)6Z(s) K|n(s)|,Z^0 A 



s,p z ,k,l * s ,p z ,k,l \ > 



s,p z ,n(s),l 



/ d P* Yl 

J s=±l 



n{s)eZ(s) l<\n(s)\,l^0 k 



3.5 The case e = — 1. 

We let h + = h+(s, I) and h- = h-(s, I) denote the differential operation h with the respective 
e = 1 and e = — 1. We then have 

h_(s, I) = io 2 d p - (7p/2 + p -1 ^) a 1 - sMa 3 = 
= ia 2 [ia 2 d p + ( 7 p/2 + p^x)^ 1 + sMa 3 ] (ia 2 ) + = ia 2 h + (-s, I) (ia 2 ) + . 

It follows that solutions = F_(s,Z, E-(s); p) of the equation (h_ — E-(s)) F_ = are 
bijectively related to solutions F + = F + (s, I, , E + (s); p) of the equation (h + — E + (s)) F + = 
by 

F_(s, /, p) = kt 2 F + (- S , Z, p), = 
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